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Foreword 


The  Betasonde  was  originally  developed  by  Pananaetrics,  Inc.  under 
the  sponsorship  of  the  Division  of  Isotopes  Development  of  the  U.  S.  A.  E.  C. 
(Contract  Monitor:  Mr.  John  C.  Dempsey)  and  the  Atmospheric  Sciences 
Branch  of  O.  N.  R.  (Scientific  Officer:  Mr.  James  Hughes).  Several 
launches  of  these  early  versions  of  the  Betasonde  were  carried  out  at 
White  Sands  Missile  Range,  under  the  supervision  of  Mr.  Harold  N. 
Ballard  of  the  Atmospheric  Sciences  Laboratory. 


The  advanced  model  Betasonde  work  reported  here  is  supported  by 
the  Army  Research  Office  and  is  monitored  by  Dr.  Arthur  V.  Dodd  of  the 
Geosciences  Division.  The  objective  of  the  work  is  to  develop  a flight 
model  of  the  Betasonde  for  integration  into  an  Areas  rocket  payload  and 
application  in  a middle  atmosphere  flight  program  at  White  Sands  Missile 
Range. 
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\k  ‘ 1.  Inti'oduction 

V 

The  Scope-of-Work  for  the  subject  contract  is  as  follows: 

The  Contractor,  as  an  independent  Contractor  and  not  as  an  agent  of 
the  Government,  shall  provide  the  necessary  management,  facilities,  ser- 
vices and  materials,  as  required,  to  accomplish  the  research  project  en- 
titled "ATMOSPHERIC  DENSITY  MEASUREMENT  IN  THE  MIDDLE  ATMOS- 
PHERE" in  accordance  with  the  Contractor's  proposal  dated  November  1975 
and  revised  during  April  1976  which  is  incorporated  as  part  of  this  contract 
by  reference.  The  Contractor  shall  use  the  types  of  personnel  listed  in  Ex- 
hibit A,  attached,  at  approximately  the  level  of  effort  stated.  The  research 
to  be  conducted  will  include  but  will  not  necessarily  be  limited  to  the  design 
and  fabrication  of  an  instrument  utilizing  beta  radiation  to  measure  atmos- 
pheric density  to  an  elevation  of  100  kilometers.  The  first  phase  of  the  re- 
search involves  investigation  of  satisfactory  beta  radiation  sources  and  de- 
tection. Space  simulation  chambers  will  be  used  to  measure  and  calculate 
the  density  response  to  beta  radiation.  Based  on  the  results,  a density 
sonde  will  be  designed  and  fabricated  in  the  second  phase  of  the  work. 

Thus,  the  work  under  this  contract  is  divided  into  two  phases: 

I Phase  I - Research  and  Development 

' ^ Phase  II  - Design  and  Fabrication 

If  use  is  to  be  made  of  the  Betasonde  fabricated  under  Phase  II,  a third 
phase  for  flight  and  Data  Analysis  will  have  to  be  funded. 

The  objective  of  the  present  report  is  to  summarize  the  Phase  I 
effort  that  has  taken  place  during  the  first  year  of  the  contract. 


2.  System  Design 

For  Phase  I,  which  is  research  and  development  work  to  determine 
optimum  performance  at  altitudes  up  to  89  km,  the  system  is  designed  for 
greatest  flexibility.  It  consists  of  a basic  betasonde  with  a charge  preamp- 
lifier and  of  modular  electronic  units  which  are  separate  from  the  sonde. 

Since  the  modules  have  adjustable  transfer  functions,  the  system  pefform- 
ance  can  be  optimized  very  efficiently.  The  source-detector  distance  of  the 
sonde  is  also  adjustable  to  provide  further  flexibility. 

The  basic  sonde  configuration  is  shown  in  Fig.  2.  1.  A tripod-supported 
boom  carries  the  annular  beta  source  which  irradiates  the  air  volume  in  the 
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region  between  source  and  detector,  outside  the  shield.  The  forward-scat- 
tered beta  rays  are  detected  by  a surface  barrier  type  semiconductor  detec- 
tor located  at  the  base  of  the  tripod.  A shield  is  placed  between  the  source 
and  detector  to  prevent  direct  passage  of  beta  particles  from  the  source  to 
the  detector,  A plastic  disc  is  placed  on  the  side  of  the  shield  that  faces  the 
source,  in  order  to  reduce  bremsstrahlung  generation.  The  distances  of  the 
source  to  the  detectors  and  to  the  shield  are  independently  adjustable  to  allow 
optimization  of  the  system  response  characteristic.  Included  in  this  develop- 
mental sonde  is  a charge  sensitive  preamplifier  which  exhibits  an  extremely 
low  noise  haracteristic.  Typically  the  noise  is  about  14  keV  FWHM  for  a 
detector  capacitance  of  200  pf  (approximate  capacitance  of  detector  to  be 
used). 

The  ^etector  is  a ruggedized  surface  barrier  type  with  an  active  area 
of  450  mm  , sufficient  to  give  a reasonable  count  rate  up  to  an  altitude  of 
80  km.  Since  a standard  detector  inherently  responds  to  light  in  the  visible 
spectrum»it  is  necessary  for  this  application  to  have  a specially  processed 
detector  featuring  a light  tight  front  surface.  According  to  the  manufacturer 
(Ortec,  Inc.)  a light  tight  front  surface  can  be  obtained  with  120pg/cm^  of 
aluminum,  which  is  vapor  deposited  directly  on  the  detector  material.  The 
latter  is  p-type  silicon  (rather  than  the  standard  n-type)  in  order  for  the 
aluminum  to  become  the  rectifying  contact  of  the  surface  barrier  diode,  and 
thus  the  front  of  the  detector.  Other  features  of  the  detector  are  low  noise 
(18  keV  max  for  beta  particles)  and  a depletion  depth  of  300p.m  at  a bias  of 
lOOV.  The  relatively  low  bias  voltage  simplifies  the  design  and  assures  low 
power  dissipation  of  the  bias  power  supply. 

The  original  concept  of  the  high-altitude  Betasonde  utilized  an  annular 
surface  barrier  detector  with  the  source  boom  mounted  through  the  center 
of  that  detector.  Although  this  arrangement  has  the  advantage  of  simpler 
mechanical  construction,  it  has  been  found  that  the  cost  of  a ruggedized, 
light  weight,  annular  detector  outweighs  these  advantages.  The  choice  of 
the  above  described  non-annular  detector  is  based  on  its  cost  effectiveness 
and  on  the  manufacturer’s  experience  with  the  construction  of  such  a special 
detector. 

The  output  signal  from  the  sonde's  charge  sensitive  preamplifier  is 
processed  by  modular  electronics  external  to  the  vacuum  chamber.  A 
block  diagram  of  the  laboratory  :.ystem  is  shown  in  Fig.  2.2,  The  voltage 
amplifier  increases  the  pulse  amplitude  of  the  sonde's  output  signal  to  con- 
venient values  by  means  of  its  adjustable  gain  control.  The  trigger  level  of 
the  threshold  detector  is  set  at  a voltage  equivalent  to  about  2,5  times  the 
energy  noise  level  of  the  detector.  As  stated  above,  the  maximum  noise 


Figure  2.2  Block  Diagram  of  Betascnde  II  Laboratory  System 
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output  from  the  detector  is  18  keV  which  results  in  an  energy  threshold 
level  of  45  keV.  The  sources  suitable  for  investigation  during  this  phase 
are:  Pm-147,  Sr-Y-90,  Ni-63  and  Tl-204.  All  have  a high  enough  maxi- 
mum energy  that  each  spectrum  contains  a useful  fraction  of  betas  above 
the  45  keV  threshold.  The  uniform  pulses  from  the  threshold  detector 
are  counted  by  the  scaler.  Thus,  count  rate  data  are  obtained  as  a function 
of  vacuum  chamber  pressure  or  by  simple  computation  as  a function  of 
density. 

For  purposes  of  detector  evaluation  a multichannel  analyzer  (MCA) 
is  included  in  the  system.  The  MCA  is  used  to  determine  the  various 
source  spectra.  Available  laboratory  type  power  supplies  are  used  to 
power  the  system. 
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3.  Optimum  Analytical  Fit  to  Altitude-Density  Profile 
3.1  General  Objective 

As  is  well  known,  density  p{z)  is  approximately  an  exponentially 
decreasing  function  of  altitude  z: 


p(z)  = 


(3.1) 


where  H is  the  scale  height  and  Pg  the  value  of  p at  z=0.  In  fact,  such  an 
approximation  is  often  made  use  of  (Ref.  3.1,  p 67)  to  provide  analytical 
representations  of  models,  for  use  in  applications  where  such  an  approach 
is  more  convenient  than  is  a long  tabulation,  for  example  where  a compu- 
tation is  to  be  made  that  requires  a continuous  repi-esentation  versus  alti- 
tude. The  result  is  generally  provided  in  the  form  of  the  coefficients  in  a 
power  series  expansion  for  fnp(z).  We  believe  that  a similar  approach 
could  also  be  useful  for  representing  the  data  from  a single  determination 
of  a p(z)  profile,  in  particular,  from  a rocket  launch  such  as  that  of  Beta- 
sonde  II,  the  instrument  being  developed  under  the  present  contract  for  air 
density  measurement  up  to  at  least  80  km  by  beta  ray  forward  scattering. 

The  Low-Background  Betasonde  (LBB)  is  an  earlier  version  Areas 
borne  instrument  designed,  principally,  for  measurement  of  p(z)  between 
about  25  and  65  km.  The  direct  output  of  the  data  reduction  procedure  is 
a tabulation  of  individual  densities  P^{z)  and  associated  statistical  uncer- 
tainties, Apj(z),  versus  z^.  When  making  use  of  data  of  this  type  it  is 
normally  not  possible  simply  to  interpolate  directly  between  values  in  the 
tabulation.  The  reason  is  that  the  points,  when  plotted,  will  generally 
show  an  amount  of  scatter  consistent  with  the  uncertainties.  Thus,  in 
order  to  use  such  data  it  is  necessary  to  draw  a best-fit  curve  of  some 
type  on  the  plot  (usually  simply  by  eye).  Values  are  then  read  off  the 
curve  and  assumed  to  have  an  uncertainty  consistent  with  that  of  the  sta- 
tistical Ap£(z),  and  any  additional  systematic  errors,  in  that  altitude  re- 
gion. Clearly,  then,  determination  of  an  optimum  analytical  fit  for  a tab- 
ulation would  serve  two  purposes.  It  would  provide  a "formal”  method  of 
fitting  the  data,  and  the  resulting  analytical  fit  would  provide  a convenient 
representation  for  computational  purposes. 

We  will  use  the  results  of  a launch  of  the  LBB  at  White  Sands  Mis- 
sile Range  at  17:00:00  UT  on  April  23,  1976  (Ref.  3.2),  to  show  the  type 
analytical  fit  that  can  be  obtained.  First  the  general  mathematical  ap- 
proach is  presented,  and  this  is  followed  by  the  density  tabulation  and  the 
application  of  the  approach  to  the  tabulation  to  determine  the  optimum  fit. 
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3.2  Mathematical  Approach 


a«piassB!^sra^g3»gqsr 


The  deviation  of  (3.1)  from  an  exact  exponential  can  be  taken  into 
account  by  allowing  the  scale  height  H to  be  a function  of  altitude.  Over 
the  altitude  region  up  to  about  100  km,  H is  not  strongly  dependent  on  z 
and  should  be  representable  as  a simple  function  (this  is  essentially  what 
is  done  in  the  standard  approach  to  model  fitting).  By  doing  this  the  re- 
sult H(z)  is  a physically  meaningful  parameter,  although  for  a given  anal- 
ytical fit  the  dependence  on  z can  be  expected  to  agree  with  model  results 
in  only  an  approximate  manner.  For  purposes  of  determining  a fit  over  a 
particular  altitude  region  it  is  then  convenient,  and  appropriate,  to  use 


p . p^e 


where  z^  is  the  (specified)  minimum  altitude  of  interest,  is  the  value 
of  p at  z = Zq  and  H(z)  is  an  "average"  value  of  H(z),  as  determined  by  the 
fitting  procedure.  Thus,  once  p©  and  H(z)  are  found  from  the  procedure, 
the  analytical  result  (3.2)  represents  the  entire  array  of  measured  density 
profile  data,  which  are  in  the  form  pj  and  statistical  uncertainty  Apj  for  all 
N altitudes  z^.  The  procedure  is  used  to  minimize  the  sum  of  the  mean 
squared  fractional  deviations  Dj  of  all  p values,  calculated  from  (3.2),  from 
the  measured  results  pj: 

Df  = (p  -Pj)/p  = 1 -p^/p  (3.3) 

Because  H(z)  occurs  in  the  exponential,  an  exact  minimi.'’ation  of  this  type 
would  require  solution  of  an  array  of  non-linear  equations.  As  we  shall 
find,  the  maximum  absolute  values  of  Dj  arc  about  0.1,  in  which  case  a 
very  good  approximation  is 

<n(p/pi)~  i - p/p  (3-4) 

For  |Dj|^|s.l,  the  difference  between  the  left  and  right  hand  sides  of  this 
equation  never  exceeds  5%. 


Thus,  we  have  used 


D.  = ln(p/p.) 


S' 


This  expression,  in  itself,  is  completely  valid  for  deriving  a set  of  values 
of  Pq.  and  H(z)  that  provides  a best  fit  to  the  experimental  data;  additionally, 
it  gives  a result  that  is  very  close  to  that  which  would  be  obtained  from  the 
exact  fit,  weighting  large  («.!)  deviations  slightly  less  {a5%)  than  does  the 
deviation  (3,3). 
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In  ordei*  to  determine  how  well  the  fit  has  been  made,  once  the  pro- 
cedure is  complete,  it  is  useful  to  compare  two  quantities.  The  first  is 
the  rms  uncertainty  the  N values  of  fractional  uncertainty  Aj^  in  the 

individual  measurement, 


2 V 2/m 

(3.6) 

where 

A;  = APj/pj 

(3.7) 

and  A pi 
is 

is  the  statistical  uncertainty  (Ref,  3.2,  p 18). 

j 

The  second  quantity 

a/  = E(Pi/p(^i)  - 

(3.8) 

which  is  the  rms  deviation  of  the  p(zj)  values,  calculated  by  use  of  (3,2), 
from  the  measured  values  p. . If  all  measurements  were  made  at  only  one 
value  of  p,  then  it  would  be  expected  that  Aj^KA^..  This  is  because,  for  a 
given  actual  p,  the  random  variations  in  the  measured  p^  for  the  Betasonde 
are  due  to  the  statistical  nature  of  the  nuclear  decay  process  that  leads  to 
beta  emission.  The  values  of  p^  are  distributed  in  a Poisson-like  manner, 
for  which  the  fractional  standard  deviation  of  the  mean  of  a series  of  N 
readings  would  be  given  by  (3.6).  The  Pearson  "chi-squared"  test  can 
readily  be  applied  in  such  circumstances  to  determine  how  well  the  mea- 
surements fit  the  distribution  assumed.  Here,  however,  we  are  dealing 
with  the  comparison  of  a best  analytical  fit  (rather  than  the  mean  of  a series 
of  measurements)  to  the  individual  measurements,  each  of  which  is  expected 
to  be  slightly  different  because  z is  changing.  Additionally,  some  atmos- 
pheric variability  will  exist  that  cannot  be  fitted  by  a simple  function  such 
as  (3,2)  without  an  unrealistically  large  amount  of  variation  in  H(z),  and 
this  will  in  itself  cause  A^,  to  exceed  A^n*  Clearly,  in  this  case  we  expect 
Ac>Am>  conclude  that  the  fit  is  satisfactory  if  the  difference  is 

acceptably  small. 

Thus,  the  procedure  here  is  to  find  the  simplest  possible  represen- 
tation of  H(z)  that  will  produce  a value  of  A^  that  is  acceptably  close  to  Aj^, 
A fit  in  which  Ac«Ani  would  not  represent  a true  best  fit;  this  could  be  ob- 
tained by  using  r sufficiently  complicated  function  for  H(z)  and  would  cor- 
respond to  following  the  random  deviations  in  the  measurements,  rather 
than  the  true  average  variation  of  p(z).  The  question  of  exactly  how  closely 
Ac  must  approach  Aj^i  for  the  analytical  fit  to  represent  the  experimental 
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tabulation  acceptably  is  not  considered  further  here.  Rather,  we  shall 
.show  that  fits  can  be  obtained  in  which  the  difference  is  small  and  reduces 
in  the  manner  expected  as  the  number  of  Aj  increases. 

I'or  use  in  the  least  square  process,  the  quantity  from  (3.5)  is 
conveniently  written 


D.  = fn(l/p.)  - fn(l/p) 


We  will  represent  H(z)  by  an  inverse  power  series  with  coefficients  Aj,  so 
that  from  (3.2)  we  have. 


fn(l/p)  = fn(l/p^)  + (z  - z^)/H(z) 


(3.10) 


wher  e 


= Ea.(z- 
j=l 


in(l/p^)  = Aj 


(3.11) 


(3.12) 


(z-  z^)/H(z) 


m 

= Da  (z-  z^)^ 

j=2 


(3.13) 


and  is  the  number  of  coefficients  A.  chosen  to  be  used.  Clearly,  once 
all  A:  in  the  single  power  series  (3.11)^are  known,  both  p^^  and  H(z)  can  be 
found  from  (3.12)  and  (3.13): 


-A  1 
e ^ 


(3.14) 


H(z)  = 1/  Ea.(z  - z )J 

/ j=2  •* 


(3.15) 


Hence,  a fit  must  be  made  to  fn(l/p)  to  determine  the  Aj  from  (3.11)  for 
use  in  (3.14)  and  (3.15),  which  then  defines  p(z)  from  (3.2). 


The  approach  used  is  a standard  least  square  procedure.  The  total 
squared  deviation  for  all  N measured  values  of  p^,  and  associated  uncer- 
tainties Apj^,  at  altitude  is  given  by  use  of  (3.9)  and  (3.10)  as 


N 


f 

i=l  ^ 

V 

1 

f 

(3.16) 

f 

1} 

where 

h.  = z.  - z 

1 1 o 

(3.17) 

?■' 

and 

f.  = fn(l/p.). 

(3.18) 

The  weights  vv.  are  taken  as  inversely  proportional  to  the  squared  fractional 
uncertainties  Aj,  ~ 


w 


, = 1/(WA.‘^) 


(3.19) 


where  Aj  is  given  by  (3.7)  and  the  normalization  constant  W is  chosen,  for 
convenience,  so  that 


= 1. 

which  defines  W from  (3.19)  as 

w = E{i/a.2). 


(3.20) 


(3.21) 


The  weighting  (3.19)  thus  gives  most  importance  to  those  measurements  for 
which  the  fractional  uncertainties  are  small;  but  none  of  the  Aj,  may  be  zero, 
of  course,  since  this  would  imply  th^-t  that  particular  is  known  exactly. 


The  least  square  procedure  then  proceeds  by  requiring  that  for  each 
particular  one  of  the  Aj,  say  Aj^, 


= O;  k = 1 to 


(3.22) 
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which,  by  use  of  (3.16),  yields  linear  equations  for  the  A,  of  the  form: 


J 


J 


m 


\ = E k = 1 to  j 
j=l 


m 


(3.23) 


where. 


N 

V,  = 'Tw.f.h. 
k -^111 
i=l 


k-1 


(3.24) 


and 


N 

B,  . = Y1 

“J  jt,  ‘ ' 


k+j-2 


(3.25) 


Thus,  given  a set  of  N values  of  p^,  A and  z^,  the  quantities  f^,  h^  and  Wj 
are  found  from  (3.18),  (3.17)  and  (3.19),  the  number  of  coefficients  Aj 
to  be  used  for  the  fit  is  chosen,  and  the  values  of  Vj^  and  are  determined 
from  (3.24)  and  (3.26)  for  k and  j = 1 to  These  values  are  then  used  in 
the  ij^  equations  (3.23),  and  the  Aj  are  determined  by  solution  of  that  set  of 
linear  equations. 

Model  fitting  often  requires  use  of  10  to  14  coefficients  in  a similar 
procedure  to  obtain  what  is  considered  to  be  an  adequate  fit.  Here  we  will 
find  that  for  the  particular  experimental  data  used,  4 or  5 are  sufficient  to 
provide  a fit  whose  accuracy  is  consistent  with  the  accuracy  of  the  data. 
Thus,  after  the  fitting  procedure  is  completed  for  each  particular  value  of 
the  value  of  A^.  is  found  from  (3.8)  and  compared  with  A^^i  from  (3.6)  to 
determine  the  adequacy  of  the  fit. 

3.  3 Flight  Data  and  Density  Tabulation 

Details  of  the  Low  Background  Betasonde  (LBB)  flight  on  April  23, 
1976  are  contained  in  Ref.  3.2  and  presented  here  only  to  the  extent  neces- 
sary for  application  of  the  above  procedure. 

The  Areas  was  launched  to  an  apogee  altitude  near  65  km,  with  the 
parachute  borne  LBB  being  ejected  just  before  apogee.  The  parachute  has 
little  effect  above  50  km.  The  descent  velocity  changes  from  about  ,2  km/ 
sec  to  .02  km/sec  between  60  and  30  km.  However,  the  sonde  does  spend 
nea?*ly  20  seconds  within  1 km  of  appgee,  thus  providing  much  data  in  this 
region  before  beginning  its  rapid  descent.  This  is  an  important  fact. 
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because  it  allows  the  density  profile  to  be  "tied  down"  at  high  alf.^ude. 
Clearly,  as  the  sonde  descent  rate  slows,  it  obtains  progressively  more 
data  in  each  altitude  region.  Hence,  the  entire  density  profile  for  the  low 
altitude  region  is  very  accurately  defined.  The  procedure  described  above 
is  designed  such  that  it  provides  a single  analytical  fit  that  automatically 
takes  all  of  these  factors  into  account. 


The  Betasonde  pulses,  each  of  which  (except  for  cosmic  ray  back- 
groxmd)  corresponds  to  detection  of  an  air- scattered  beta  particle,  are 
telemetered  and  recorded  on  magnetic  tape  in  real-time.  These  pulses 
are  then  summed  into  approximately  equal  time  intervals  tj^  - depending 
upon  telemetry  noise  and  time  used  for  transmitting  temperature  for  di- 
agnostic purposes.  In  Table  3.1  these  sums  are  labelled  CTS  at  each 
altitude  2(km)  for  the  time  interval  TM(SECS).  These  are  the  basic  flight 
data  from  which  the  density  measurement  at  each  altitude  is  determined. 


It  is  seen  that  the  counts  generally  increase  as  altitude  decreases, 
because  the  number  of  beta  particles  scattered  increases  as  the  density 
increases  in  this  density  region.  In  addition  lo  the  density  itself,  there 
are  several  other  factors  that  must  be  taken  into  account  in  determining 
density  from  this  count  rate; 


1)  circuit  dead-time,  caused  by  dead  time  in  geiger  counters  after 
counting  a pulse  (KlSfisec)  and  a dead-time  circuit  introduced 
purposely  to  keep  the  count  rate  into  the  telemetry  transmitter 
from  exceeding  its  capabilities, 

2)  cosmic  ray  effect,  a small,  but  non-neglig\ble,  altitude  depen- 
dent effect, 

and  3)  time  elapsed  between  calibration  and  flight. 


Each  of  these  factors  has  been  taken  into  account  as  necessary  in  obtain- 
ing the  calibration  data  of  Table  3.2  and  the  flight  density  tabulation  in 
Table  3.3.  In  the  latter  table  the  count  data  from  Table  3.1  are  repeated 
as  the  first  three  columns,  with  the  observed  counts  (OBSCTS)  in  the  third. 
This  leads  to  an  observed  count  rate  (not  shown)  which  is  corrected  ^or 
dead  time,  by  use  of  the  circuit  dead  time  of  185;isec,  to  obtain  the, correc- 
ted count  rate  (CORCTR).  The  cosmic  ray  count  rate  (COSCTR)  is  then 
subtracted  to  obtain  the  density-dependent  count  rate  (DENCTR).  The 
DENSITY  (g/m^)p^  for  the  altitude  is  then  determined  by  use  of  thp 
calibration  results  of  Table  3.2,  and  the  absolute  uncertainty,  UNC(g/tn^) 
from  the  density,  the  counts,  and  the  calibration  data.  The  fractional  un- 
certainty is  given  by  the  ratio  UNC/DENSITY. 
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Table  3.1  Low  Background  Betasonde  Count  Data  for  April  23,  1976  Flight 


Z(KM> 

TM<SCCS> 

CIS 

64.0 

5.5 

80. 

64.3 

5.2 

89. 

64.2 

5.4 

85. 

63.5 

5.2 

83. 

62.9 

5.1 

102. 

62.1 

5.1 

99. 

61 .0 

5.2 

1 30. 

60.0 

5.2 

1 33. 

59.1 

5.1 

156. 

57.9 

5.1 

155. 

55.9 

5.1 

219. 

54.9 

5.1 

240* 

53.9 

5.1 

259. 

53.0 

5.1 

279. 

52.0 

5.0 

355. 

51  .2 

5.2 

376. 

50.5 

5.0 

376. 

50.0 

4.1 

31  1 . 

49.5 

4.0 

361  . 

48.6 

5.0 

543. 

48.1 

5.0 

502. 

47.7 

5.1 

565* 

47.2 

4.0 

467. 

46.8 

4.0 

485. 

46.  4 

4.1 

495. 

45.8 

2*0 

307. 

45.3 

4.0 

650. 

45.0 

4*0 

626. 

44.2 

4*0 

683* 

43.9 

4.0 

748. 

43.4 

5.0 

1047. 

43.1 

5.1 

989. 

42.8 

5.0 

1093. 

42.5 

5.0 

1124. 

42.2 

5.1 

1 197. 

41  .9 

4.0 

942. 

41.7 

4.0 

1024. 

41.1 

5.0 

1408. 

40*8 

S.O 

1 443* 

40.5 

5.0 

1510. 

40.2 

5.1 

1538. 

40.0 

5.0 

1610. 

39.7 

5.0 

1652. 

39.5 

5.0 

1742. 

39.3 

5.0 

1803. 

39.0 

4.0 

1507. 

38.5 

4.  1 

1574. 

38.3 

5.0 

2007. 

38.1 

5.0 

208  6. 

37.9 

5.0 

2228. 

37.6 

5.0 

2241  * 

37.4 

5.0 

2335. 

37.2 

5.0 

2385. 

37.0 

5.0 

2400. 

36*8 

5.0 

2545. 

36.2 

S.O 

2808. 

36.1 

5.0 

2964. 

35.9 

S.O 

2975. 

35.7 

5.0 

2890. 

35.5 

5.0 

3102. 

35.4 

S.O 

3098. 

35.2 

5.0 

3166. 

35.1 

5.0 

3344. 

34.9 

4.0 

2655. 

34.5 

5.1 

3547. 

34.4 

5.0 

3628. 

34.3 

5.0 

3530. 

34.1 

5.0 

3444. 

34.0 

5.0 

3666. 

33.9 

5.0 

3809. 

33.7 

S.O 

3829. 

33.6 

S.O 

4038. 

33.5 

4.0 

3283* 

33.1 

5.0 

4097. 

33.0 

5.1 

4227. 

32.9 

5.0 

4359. 

32.8 

5.0 

4363. 

32  .'6 

5.0 

4580. 

32.4 

5.0 

4735. 

32.2 

S.O 

4755. 

32.1 

5.0 

4895. 

32.0 

4.0 

3944. 

31 .6 

5.0 

5055. 

31.5 

5.0 

5353. 

31.3 

5.1 

5316. 

31  .2 

5*0 

5271  . 

31  .0 

5.0 

5406. 

30.8 

5.0 

5701  . 

30.7 

5.0 

5707. 

30.5 

5.0 

5901  . 

30.4 

4.0 

4735. 

30.0 

5.0 

6338. 

29.9 

5.0 

6418. 

89.7 

5.0 

6332. 

29.6 

5.0 

6350. 

29.5 

5.0 

6462. 

29.3 

5.0 

6764. 

29.1 

5.0 

6708* 

29.0 

5.0 

6757. 

28.9 

4.0 

5757. 

28.6 

5.0 

6895. 

28.5 

5.1 

7383,. 

28.3 

5.0 

7335. 

28.2 

5.0 

7454. 

28.1 

S.O 

7724. 

27.9 

5.0 

7764. 

27.7 

5.0 

7901. 

27.6 

5.0 

7778. 

27.4 

4.0 

6502* 

27.2 

5.0 

8424* 

27.1 

4.0 

6541  . 

27.0 

5.0 

8401  . 

26.9 

5.0 

8239. 

26.7 

5.0 

8719. 

26*6 

5.0 

8675. 

26.5 

5.1 

8509. 

26.3 

5*0 

8914. 

26.2 

S.O 

8892. 
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Table  3.2  Low  Background  Betasonde  Density  Calib- 
ration Data  Recorded  April  7,  1976. 


01 

00 

0000 

02 

02 

01  74 

03 

04 

0352 

04 

06 

0524 

05 

08 

0700 

06 

10 

08  74 

07 

12 

1045 

08 

14 

1220 

09 

16 

1400 

10 

18 

1575 

1 1 

20 

1750 

12 

22 

1900 

1 3 

24 

2020 

14 

26 

21  30 

15 

28 

2245 

16 

30 

2345 

17 

32 

2440 

18 

34 

2530 

19 

36 

2620 

20 

38 

2705 

21 

40 

2790 

22 

42 

2880 

Table  3.3  Low  Background  Betasonde  Measured  Density 
Tabulation  for  April  23,  1976  Flight 

M EAN-L1FE« DAYS” 1 378. 

TIME  FROM  CAL# DAYS*  16. 

CIR.  DEAD  TIME#SEC>. 000185 

UNITS* ALTaKMI TIME-INT»SECl ALL  CTR*CPS*  DEN# UNC=G/M**  3 


ALT 

TIME-INT 

OBSCTS 

CORCTR 

COSCTR 

DENCTR 

DENSITY 

UNC  FRAC-UNC 

64.0 

5.5 

80. 

14.58 

0.30 

1 4.45 

0.168 

0.019 

0.112 

64.3 

5.2 

89. 

17.17 

0.30 

17.07 

0.198 

0.021 

0.106 

64.2 

5.4 

85. 

15.79 

p.30 

15.67 

0*182 

0.020 

0.108 

63.5 

5.2 

83. 

16.01 

0.30 

15.89 

0.185 

0.020 

0.110 

62.9 

5.1 

102. 

20.07 

0.30 

20.01 

0.232 

0.023 

0.099 

62.1 

5.1 

99. 

19.48 

p.30 

19.41 

0.225 

0.023 

0.100 

61  .0 

5.2 

130. 

25.12 

p.30 

25.1  1 

0.291 

0.026 

0.088 

60.0 

5*2 

1 33. 

25.70 

0.30 

25.70 

0.298 

0.026 

0.087 

59.1 

5.1 

156. 

30.76 

6.30 

30.82 

0. 357 

0.029 

0.080 

57.9 

5.1 

155. 

30.56 

0.30 

30.62 

0.355 

0.028 

0.080 

55.9 

5.1 

219. 

43..29 

6.30 

43.49 

0.504 

0.034 

0.067 

54.9 

5.1 

240. 

47.47 

0.30 

47.72 

0.553 

0.036 

0.064 

53.9 

5.1 

259. 

51  .27 

0.30 

51.56 

0.597 

0.037 

0.062 

53.0 

5.1 

279. 

55.27 

0.30 

55.61 

0.644 

0.038 

0.060 

52.0 

5.0 

355. 

71  .94 

0.30 

72.48 

0.8  39 

0.044 

0.053 

51.2 

5.2 

376. 

73.29 

0.30 

73.84 

0.854 

0.044 

0.051 

50.5 

5.0 

376. 

76.26 

0.30 

76.85 

0.889 

0.046 

0.051 

Table  3.  3 (Cont,) 


I 


fe- 

I 

£ 


i 


I 


I 


ft- 


i- 

f 

I 

y 

f 

!■ 

f- 


[' 

X 'A 


ALT 

TIME-INT 

OBSCTS 

CORCTH 

COSCTR 

DENCTR 

DENSITY 

UNO  FRAC-UNC 

SO.O 

4.1 

31 1 . 

76.93 

0.30 

77.53 

0.897 

0.051 

0.056 

49  • 5 

4.0 

361  . 

91 .78 

0.31 

92.54 

1 .069 

0.056 

0.052 

48  • 6 

5.0 

543. 

1 10.83 

0.32 

1 1 1 .80 

1 .290 

0.055 

0.043 

48*1 

5.0 

502. 

102.30 

0.32 

103.1  7 

1.191 

0.053 

0.044 

47.7 

5.1 

565  . 

113.10 

0.33 

114.09 

1 .316 

0.055 

0.042 

47.2 

4.0 

467. 

119.33 

0.33 

1 20 . 38 

1.388 

0.064 

0.046 

46*8 

4.0 

485. 

124.03 

0.34 

125.14 

1.443 

0.065 

0.045 

46.  4 

4.  1 

495. 

123.49 

0.34 

124.58 

1.437 

0.064 

0.045 

4S  . 8 

2.0 

307. 

157.99 

0.35 

159.48 

1 .835 

0.104 

0.056 

45*3 

4.0 

650. 

167.54 

0.36 

169.1  3 

1 .945 

0.075 

0.039 

45  . C 

4.0 

626. 

161 .1  7 

0.36 

162.68 

1 .871 

0.074 

0.040 

44.2 

4.0 

683. 

176.32 

0.37 

178.00 

2.044 

0.075 

0.037 

43.9 

4.0 

748. 

193.70 

0.37 

195.58 

2.239 

0.079 

0.035 

43. 4 

5.0 

1047. 

217.84 

0.38 

220.00 

2.510 

0.076 

0.030 

43.  1 

5*1 

989. 

201.14 

0.38 

203.10 

2,322 

0.072 

0.031 

42 . 8 

5.0 

1093. 

227.81 

0.39 

230.08 

2.622 

0.078 

0.030 

42.5 

5.0 

1124. 

234.55 

0.39 

236.90 

2.699 

0.079 

0.029 

42 . 2 

5.1 

1197. 

245.36 

0.40 

247.82 

2.821 

0.080 

0.028 

41  .9 

4.0 

942. 

246.23 

0.40 

248 . 70 

2.831 

0.091 

0.P32 

41.7 

4.0 

1024. 

268.73 

0.40 

271 .46 

3.086 

0.095 

0.031 

41 . 1 

5.0 

1408. 

297.08 

0.41 

300.13 

3.410 

0.091 

0.027 

40.8 

5.0 

1443. 

304.88 

0.41 

308.02 

3.499 

0.092 

0.026 

40.5 

5.0 

1510. 

319.87 

0.42 

323.19 

3.671 

0.095 

0.026 

40.2 

5.1 

1538. 

319.39 

0.42 

322.69 

3.666 

0.094 

0.026 

40.0 

5.0 

1610. 

342.40 

0.42 

345.97 

3.931 

0.098 

0.025 

39.7 

5.0 

1652. 

351 .91 

0.43 

355.59 

4.042 

0.103 

0.025 

39 . 5 

5.0 

1742. 

372.40 

0.43 

376.32 

4.285 

0.106 

0.025 

39.3 

5.0 

1803. 

386.38 

0.43 

390.45 

4.451 

0.108 

0.024 

39.0 

4.0 

1507. 

404.98 

0.43 

409.27 

4.671 

0.123 

0.026 

38 . 5 

4. 1 

1574. 

413.25 

0.44 

41  7.63 

4.768 

0.123 

0.026 

38.3 

5.0 

2007. 

433.60 

0.44 

438.21 

5.008 

0.114 

0.023 

38 . 1 

5.0 

2086. 

452.09 

0.45 

456.92 

5.225 

0.116 

0.022 

37.9 

5.0 

2228. 

485.63 

0.45 

490.85 

5.618 

0.120 

0.021 

37  • 6 

5.0 

2241. 

488  • 72 

0.45 

493.97 

5.654 

0.120 

0.021 

37.4 

5.0 

2335. 

511.16 

0.45 

516.67 

5.916 

0.123 

0.021 

37.2 

5.0 

2385. 

523.1 7 

0.46 

528.81 

6.054 

0.122 

0.020 

37.0 

5.0 

2400. 

526.78 

0.46 

532.47 

6.096 

0.123 

0.020 

36.8 

5.0 

2545. 

561 .91 

0.46 

568.01 

6.498 

0.128 

0.020 

36.2 

5.0 

2808. 

626.71 

0.47 

633.56 

7.242 

0.136 

0.019 

36.  1 

5.0 

2964. 

665.82 

o;47 

673.12 

7.693 

0.141 

0.018 

35.9 

5.0 

2975. 

668.60 

0.47 

675.93 

7.725 

0.141 

0.018 

35  . 7 

5.0 

2890. 

647.21 

0.48 

654.28 

7.478 

0.1  39 

0.019 

35 . 5 

5.0 

3102. 

700.84 

O:  48 

708.54 

8.097 

0.145 

0.018 

35.4 

5.0 

3098. 

699.82 

0 . 48 

707.51 

8.086 

0.145 

0.018 

35.2 

5.0 

3166. 

71 7.22 

0^48 

725.11 

8.286 

0.147 

0.018 

35. 1 

5.0 

3344. 

763.23 

0!48 

771 .66 

8.819 

0.153 

0.017 

34.9 

4.0 

2655. 

756.66 

6.49 

765.01 

8.743 

0.170 

0.019 

34.  5 

5.1 

3547. 

798.19 

0.49 

807.02 

9.226 

0.156 

0.017 

34.4 

5.0 

3628. 

838.10 

0.49 

847.39 

9.692 

0.163 

0.017 

34.3 

5.0 

3530. 

812.06 

0.49 

821 .05 

9.388 

0.159 

0.017 

34. 1 

5.0 

3444. 

789.39 

0.50 

798.11 

9.123 

0.156 

0.017 
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Table  3.3  (Cont.) 


ALT 

TIME-INT 

OBSCTS 

COKCTK 

COSCTR 

DENCTR 

DENSITY 

UNC  FRAC-UNC 

34.0 

5.0 

3666. 

848.26 

0.50 

857.66 

9.811 

0.164 

0.017 

33.9 

5.0 

3809. 

886.78 

0.50 

896.63 

10.267 

0.171 

0.017 

33.7 

5.0 

3829. 

892.20 

0.50 

902. 1 1 

10.332 

0.172 

0.017 

33.6 

5.0 

4038. 

949.45 

0.50 

960.04 

11.012 

0.177 

0.016 

33.5 

4.0 

3283. 

967.68 

0.50 

978.47 

1 1 .227 

0.199 

0.018 

33.1 

5.0 

4097. 

965.81 

0.51 

976.57 

1 1 .205 

0.178 

0.016 

33.0 

5.1 

4227. 

978.92 

0.51 

989.84 

1 1 . 360 

0.177 

0.016 

32.9 

5.0 

4359. 

1039.44 

0.51 

1051 .07 

12.070 

0.184 

0.015 

32.8 

5.0 

4363. 

1040.58 

0.51 

1052.22 

12.084 

0.184 

0.015 

32.6 

5.0 

4580. 

1102.90 

0.51 

1 115.26 

12.810 

0.189 

0.015 

32.4 

5.0 

4735. 

1148.15 

0.52 

1 161  .04 

1 3.332 

0.192 

0.014 

32.2 

5.0 

4755. 

1 1 54.04 

0.52 

1166.99 

1 3.400 

0.192 

0.014 

32.1 

5.0 

4895. 

1195.53 

0.52 

1208.96 

1 3.876 

0.195 

0.014 

32.0 

4.0 

3944. 

1205.98 

0.52 

1219.54 

13.995 

0.219 

0.016 

31 .6 

5.0 

5055. 

1243.60 

0.53 

1257.59 

14.413 

0.195 

0.014 

31  .5 

5.0 

5353. 

1 335.01 

0.53 

1 350.07 

15.439 

0.206 

0.01  3 

31.3 

5.1 

531  6. 

1 29 1 . 38 

0.53 

1 305.92 

14.947 

0.199 

0.013 

31  .2 

5.0 

5271. 

1 309.61 

0.53 

1324.37 

15.153 

0.203 

0.013 

31.0 

5.0 

5406. 

1351 .54 

0.53 

1 366.78 

15.627 

0.208 

0.013 

30.8 

5.0 

5701 . 

1445.00 

0.54 

1461 .34 

16.701 

0.22! 

0.013 

30.7 

5.0 

5707. 

1 446.93 

0.54 

1463.29 

16.723 

0.221 

0.'013 

30.5 

5.0 

5901. 

1509.86 

0.54 

1526.95 

17.451 

0.227 

0.013 

30.4 

4.0 

4735. 

1515.67 

0.54 

1532.83 

17.518 

0.255 

0.015 

30.0 

5.0 

6338. 

1655.92 

0.55 

1674.71 

19.096 

0.243 

0.013 

29.9 

5.0 

6418. 

1683.33 

0.55 

1 702.44 

19.421 

0.253 

0.013 

29.7 

5.0 

6332. 

1653.88 

0.55 

1672.64 

19.072 

0.243 

0.013 

29.6 

5.0 

6350. 

1660.02 

0.55 

1678.85 

19.144 

0.245 

0.013 

29.5 

5.0 

6462. 

1698.50 

0.55 

1717.78 

19.605 

0.258 

0.01  3 

29.3 

5.0 

6764. 

1804.38 

0.55 

1824.89 

20.929 

0.296 

0.014 

29.1 

5.0 

6708. 

1784.51 

0.56 

1804.79 

20.666 

0.283 

0.014 

29.0 

5.0 

6757. 

1801.89 

0.56 

1822.37 

20.896 

0.294 

0.014 

28.9 

4.0 

5757. 

1961 .53 

0.56 

1983.87 

23.378 

0.444 

0.019 

28 . 6 

5.0 

6895. 

1851 .29 

0.56 

1872.34 

21.589 

0.327 

0.015 

28.5 

5.1 

7383. 

1977.16 

0.56 

1999.68 

23.648 

0.400 

0.017 

28.3 

5.0 

7335. 

201 3.44 

0.57 

2036.  38 

24.303 

0.439 

0.018 

28.2 

5.0 

7454. 

2058.54 

0.57 

2082.01 

25. 1 38 

0.43(7 

0.017 

28.1 

5.0 

7724. 

2162.94 

0.57 

2187.63 

26.962 

0.433 

0.016 

27.9 

5.0 

7764. 

21 78.66 

0.57 

2203.53 

27.242 

0.445 

0.016 

27.7 

5.0 

7901. 

2232.98 

0.57 

2258.48 

28.263 

0.49'8 

0.018 

27.6 

5.0 

7778. 

2184.17 

0.57 

2209.10 

27.341 

0.449 

0.016 

27.4 

4.0 

6502. 

2324.53 

0.58 

2351 .09 

30.125 

0.599 

0.020 

27.2 

5.0 

8424. 

2447.73 

0.58 

2475.73 

32.794 

0.599 

0.018 

27.1 

4.0 

6541 . 

2344.52 

q.58 

2371 .31 

30.542 

0.609 

0.020 

27.0 

5.0 

8401  . 

24  38.03 

Q.58 

2465.91 

32.576 

0.598 

0.018 

26.9 

5.0 

8239. 

2370.40 

C.58 

2397.49 

31  .091 

0.558 

0.018 

26.7 

5.0 

8719. 

2574.27 

0.59 

2603.74 

35.629 

0.6q2 

0.018 

26.6 

5.0 

8675. 

2555.1 3 

C.59 

2584.  38 

35.194 

0.620 

0.018 

26.5 

5.1 

8509. 

241  3.33 

0.59 

2440.92 

32.020 

0.588 

0.018 

26.3 

5.0 

8914. 

2660.17 

0.59 

2690.64 

37. 662 

0.671 

0.018 

26.2 

5.0 

8892. 

2650.39 

0.59 

2680.74 

37.429 

0.669 

0.018 

The  last  of  the  factors  that  must  be  taken  into  account  is  important 
because  of  decay  in  intensity  of  the  radioactive  source,  in  this  case  Pm-147 
whose  mean  life  is  1378  days.  Essentially,  a calibration  curve  that  would 
be  correct  on  the  day  of  the  flight  is  derived  from  that  recorded  earlier. 

In  the  present  instance  the  time  lapse  was  only  16  days,  and  the  calibra- 
tion data  are  shown  in  Table  3.2.  There  are  22  points  at  which  the  density 
p(g/m2)  and  the  count  rate,  corrected  for  the  dead  time  factor,  was  as 
shown  in  CTS/SEC.  At  low  values  of  p the  count  rate  increases  approx- 
imately linearly,  but  as  p increases  the  count  rate  increases  less  rapidly. 
At  25  km  Pr:40  g/m^  and  at  30  km  p~20  g/m^.  Data  below  about  30  km 
were  not  analyzed  in  Ref.  3.1  because  of  the  non-linear  nature  of  the  calib- 
ration curve  in  that  region.  Here,  however,  all  of  the  data  in  Table  3.1 
are  analyzed.  This  was  done  by  making,  for  each  flight  count  rate,  a 
quadratic  fit  to  three  appropriate  calibration  count  rate  points  in  Table 
3.2.  The  measured  DENSITY  values  of  Table  3.3  were  thus  determined 
from  that  analytical  fit  to  the  density  - count  rate  results  in  Table  3.2. 

The  array  of  118  altitudes  Zj,  measured  densities  p^  and  fractional 
statistical  uncertainties  Aj  given  in  Table  3.3  are  the  input  data  for  the 
least  square  mathematical  procedure  described  in  Section  3.2. 

3.4  Results  of  Analytical  Fit  for  April  23,  1976  Low  Background 
Betasonde  Flight 

By  use  of  the  data  in  Table  3.3  for  z^,  p^  and  Aj^,  Eqs.  (3.23)  were 
solved  for  - 2-5  to  obtain  the  values  for  the  coefficients  Aj.  For  each 
set  of  solutions  the  results  for  rms  sum  of  the  weighted  deviations  between 
the  analytical  lit  and  experimental  points  was  determined  from  (3.16). 

This  is  the  quantity  that  was  minimized  to  produce  the  solution  obtained. 
Ac.  the  rms  unweighted  deviation  sum,  was  also  calculated  from  (3.8); 
both  A and  Ac  are  to  be  compared  to  Aj^.  the  rms  fractional  uncertainty, 
from  (3.6).  Results  for  all  of  these  parameters,  along  with  p^  from  (3.14), 
are  given  in  Table  3.4. 

The  rms  fractional  uncertainty,  is  independent  of  and  ha." 

the  value  0.039.  As  noted  above,  a satisfactory  fit  is  one  in  which  the 
difference  between  A^,  and  A^  is  acceptably  small. 

As  seen  from  the  table,  A decreases  rapidly  as  goes  from  2 to  3, 
but  with  j'm  as  large  as  5 the  accuracy  of  the  fit  is  marginally  better.  A^ 
appears  to  go  thru  a minimum  at  = 4.  As  observed  above,  a fit  derived 
by  minimizing  A,  the  rms  siun  of  the  deviations  weighted  inversely  with  the 
uncertainties,  is  taken  here  as  the  proper  procedure  to  follow.  This 
causes  those  points  with  most  accuracy  to  produce  the  most  effect  on  the 
fit.  However,  such  a fit  should  not  be  allowed  to  produce  a large  value  of 
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n 

k 


Aj,,  the  rms  sum  of  the  unweighted  deviations.  In  vi  ;w  C'f  the  that 
there  is  undoubtedly  some  short-term  atmospheric  variability  .n  density, 
it  is  unlikely  that  a fit  in  which  approaches  much  ti.  closely 
than  in  the  present  case  would  be  expected.  Thus,  in  the  present  instance 
at  least,  a criterion  for  selecting  the  most  satisfactory  of  the  fits  obtained 
by  minimizing  A,  is  that  particular  one  for  which  A^.  is  a minimum;  = 4 
here.  Whether  such  a dependence  on  will  always  occur  is  not  presently 
known.  When  the  results  of  additional  flights  are  available,  further  inves- 
tigation of  the  method  of  selecting  the  optimum  fit  shovild  be  made. 

The  quantity  H(z)  given  by  (3.15)  is  cctually  to  be  used  in  (3.2)  to  pro- 
vide the  optimum  analytical  fit  to  the  data.  However,  for  purposes  of  com- 
parison with  models  it  is  useful  to  have  the  value  of  H(z).  This  can  be  de- 
termined as  follows.  Note  from  (3.1)  that 


H(z)  = 


But,  from  (3.2)  for  the  fit  used  here  this  gives 


Use  of  (3.15)  then  yields 

H(z) 


/ m . 

yT.  A (j-l)(z-Zo)^“ 

M=2 


(3.26) 


(3.27) 


(3.28) 


Consequently,  once  the  Aj  have  been  determined  it  is  possible  to  calc\ilate 
the  value  of  H(z)  as  well  as  H(z). 

For  the  fit  with  = 4,  results  are  given  in  Table  3.6  and  Figures 
3.1  and  3.2.  The  first  three  columns  of  the  table  give  z-,  pj  and  Apj  (ALT, 
DENSITY  and  UNC  in  Table  3.3).  The  next  three  give  p(z),  H(z),  and  H(z) 
as  calculated  from  (3.2),  (3.15),  and  (3.28),  respectively.  The  next  column 
gives  the  fractional  difference  between  measured  and  calculated  p(z). 

The  final  column,  labeled  RMS  FRA.C,  is  a running  rms  average  of  the 
fractional  difference  for  that  altitude  with  that  of  the  four  previous  (higher 
altitude  points).  Hence,  it  is  assigned  the  value  zero  for  the  four  highest 
altitude  points. 
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Table  3.5  Altitude  Dependence  Results  of 
Analytical  Fit  for  LBB  Flight  on 
April  23,  1976 


j 


m 


= 4 


results  of  fit 

J A(J> 

1 -0.35753E+01 

2 0.17017E-»00 

3 -0.76747E-03 

4 -0.18776E-05 

Z0-  26.2  KM 
RHO0>  35.707  G/M4>4  3 


RMS  FRAC-RMS  AV6  OF  PREV  5 DIFF  FRACS 


ALT 

RHO-MEAS 

UNGER 

RHO-CALC 

HBAR<Z> 

H(Z> 

DIFF 

RMS 

KM 

G/M«4t3 

G/M4»l>3 

G/M**3 

KM 

KM 

FRAG 

FRAG 

64.0 

0.168 

0.019 

0.190 

7.221 

9.606 

-0.117 

0.000 

64.3 

0.198 

0.021 

0.184 

7.236 

9.661 

0.073 

0.000 

64.2 

0.182 

0.020 

0.186 

7.23J 

9.642 

-0.024 

0.000 

63. S 

0.185 

0.020 

0.201 

7.198 

9.516 

-0.077 

0.000 

62.9 

0.232 

0.023 

0.214 

7.170 

9.412 

0.085 

0.081 

62.1 

0.225 

0.023 

0.233 

7.133 

9.276 

-0.033 

0.064 

61.0 

0.291 

0.026 

0.262 

7.083 

9.096 

0.109 

0.073 

60.0 

0.298 

0.026 

0.293 

7.038 

8.940 

0.017 

0.073 

59.1 

0.357 

0.029 

0.324 

6.998 

8.805 

0.100 

0.078 

57.9 

0.355 

0.028 

0.372 

6.947 

8.632  - 

-0.046 

0.071 

55.9 

0.504 

0.034 

0.471 

6.862 

8.360 

0.070 

0.077 

54.9 

0.553 

0.036 

0.532 

6.821 

8.232 

0.040 

0.062 

53.9 

0.597 

0.037 

0.60! 

6.781 

8.108 

-0.006 

0.061 

53.0 

0.644 

0.038 

0.672 

6.745 

8.001 

-0.041 

0.046 

52.0 

0.839 

0.044 

0.762 

6.706 

7.885 

0.101 

0.061 

51.2 

0.854 

0.044 

0.844 

6.675 

7.796 

0.012 

0.052 

50.5 

0.889 

0.046 

0.923 

6 . 648 

7.719 

-0.037 

0.05R 

50.0 

0.897 

0.051 

0.985 

6.629 

7.666 

-0.090 

0.066 

49.5 

1.069 

0.056 

1 .052 

6.61 1 

7.613 

0.016 

0.063 

48.6 

1 .290 

0.055 

1.185 

6.577 

7.521 

0.089 

0.060 

48.1 

1.191 

0.053 

1.267 

6.559 

7.471 

-0.060 

0.06^ 

47.7 

1 .316 

0.055 

1.337 

6.544 

7.431 

-0.015 

0.063 

47.2 

1.388 

0.064 

1.430 

6.526 

7.383 

-0.029 

0.051 

46.8 

1 .443 

0.065 

1 .510 

6.512 

7.344 

-0.044 

0.054 

46.4 

1 .437 

0.064 

1.594 

6.498 

7.306 

-0.099 

0.057 

45.8 

1.835 

0.104 

1.731 

6.476 

7.250 

0.060 

0.057 

45.3 

1 .945 

0.075 

1 .855 

6.459 

7.205 

0.048 

0.061 

45.0 

1 .871 

0.074 

1 .935 

6.448 

7.178 

-0.033 

0.061 

44.2 

2.044 

0.075 

2.164 

6.421 

7.106 

-0.055 

0.063 

43.9 

2.239 

0.079 

2.257 

6.410 

7.080 

-0.008 

0.045 

43.4 

2.510 

0.076 

2.423 

6.393 

7.037 

0.036 

0.040 

43.1 

2.322 

0.072 

2.529 

6.383 

7.012 

-0.082 

0.049 
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Table  3.  5 (coni’ d) 


ALT  FKO-MEAS 

UNCER  RKO-CALC 

X)$ 

a/H«43 

fl/H**3 

G/M**3 

42.8 

2.622 

0.078 

2*639 

42.5 

2.699 

0.079 

2.755 

42.2 

2.821 

0.080 

2.877 

41.9 

2.831 

0.091 

3.004 

41.7 

3.086 

0.095 

3.093 

41.1 

3.410 

0.091 

3.375 

40.8 

3.499 

0.092 

3.526 

40.5 

3.671 

0.095 

3.665 

40.2 

3.666 

0.094 

3.852 

40.0 

3.931 

0.098 

3.967 

39.7 

4.042 

0.103 

4.148 

39.5 

4.285 

0.106 

4.273 

39.3 

4.451 

0.108 

4.402 

39.0 

4.671 

0.123 

4.604 

38.5 

4.768 

0.123 

4.962 

38.3 

S.008 

0.114 

5.114 

38.1 

5.225 

0.116 

5.271 

37.9 

5.618 

0.120 

5.433 

37.6 

5.654 

0.120 

5.685 

37.4 

5.916 

0.123 

5.861 

37.2 

6.054 

0.122 

6.043 

37.0 

6.096 

0.123 

6.230 

36.8 

6.498 

0*128 

6.424 

36.2 

7.242 

0.136 

7.045 

36.1 

7.693 

0.141 

7.154 

35.9 

7.725 

0.141 

7.379 

35.7 

7.478 

0.139 

7.611 

35.5 

8.097 

9.145 

7.851 

35.4 

8.086 

0.145 

7.974 

35.2 

8.286 

0.147 

3*226 

35.1 

8.819 

0.153 

8.356 

34.9 

8.743 

0.170 

8.621 

34.5 

9.226 

0.156 

9.179 

34.4 

9.692 

0.163 

9.324 

34.3 

9.388 

0.159 

9.472 

34.1 

9.123 

0.156 

9.775 

34.0 

9.811 

0.164 

9.930 

33.9 

10.267 

0.171 

10.088 

33.7 

10.332 

0.172 

10.413 

33.6 

11.012 

0.177 

10.579 

33.5 

11.227 

0.199 

10.748 

33.1 

11.205 

0.178 

11.454 

33.0 

11.360 

0.177 

11.638 

32.9 

12.070 

0.184 

11.825 

32.8 

12.084 

0.184 

12.015 

32.6 

12.810 

0.189 

12.406 

HBAR(2> 

H(Z) 

DIFF 

RMS 

KM 

KM 

FRAC 

FRAC 

6.373 

^.986 

-0.007 

0.047 

6.363 

6.961 

-0.021 

0.041 

6.353 

6.936 

-0.019 

0.042 

6.343 

6.912 

-0.058 

0.047 

6.336 

6.895 

-0.002 

0.029 

6.316 

6.847 

0.010 

0.029 

6.307 

6.823 

-0  * 008 

0.028 

6.297 

6.800 

-0.004 

0.027 

6.287 

6.776 

-0.048 

0.022 

6.281 

6.761 

-0.009 

0.023 

6.271 

6.738 

-0.025 

0.025 

6.264 

6.722 

0.003 

0.025 

6.258 

6.707 

0.01 1 

0.025 

6.248 

6.684 

0.015 

0.015 

6.233 

6 '647 

-0.039 

0.022 

6.226 

6.632 

-0.021 

0.021 

6.220 

6.618 

-0.009 

0.022 

6v214 

6.603 

0.034 

0.026 

6.204 

6.582 

-0.006 

0.025 

6.198 

6.567 

0.009 

0.019 

6.192 

6.553 

0.002 

0.017 

6.186 

6.539 

-0.022 

0.019 

6*180 

6.525 

0.012 

0.012 

6.161 

6.483 

0.028 

0.017 

6*158 

6.476 

0.075 

0.038 

6.152 

6.1A6 

6.140 

6.137 

6.131 

6.128 

6.122 

6.110 

6.107 

6.104 

6-098 

6.095 

6.092 

6.086 

6.083 

6.080 

6.068 

6.066 

6.063 

6.060 

6.054 


6.462 

6.448 

6.435 

6.428 

6.414 

6.408 

6.394 

6.368 

6.361 

6.355 

6.341 

6.335 

6.328 

6.315 

6.309 

6.303 

6.277 

C.271 

6.264 

6.258 

6.245 


0.047 

-0.017 

0.031 

0.014 

0.007 

0.055 

0.014 

0.005 

0.039 

-0.009 

-0.067 

-0.012 

0.018 

-0.fed8 

0.041 

0.045 

-0.022 

-0.024 

0.021 

0.006 

0.033 


0.043 

0.043 

0.045 

0.043 

0.027 

0.030 

0.030 

0.087 

0.031 

0.031 

0.036 

0.035 

0.036 

0.032 

0.036 

0.029 

0.030 

0.031 

0.032 

0.026 

0.023 


\ 


Table  3.  5 (cont'd) 


A> 


^ ,»» 


ALT 

MtO'MEAS 

UNCER 

RKO^CALC 

HBARCZ) 

H».:- ) 

DIFF 

RMS 

XM 

6/M4>*  3 

9/M44>  3 

Q/M*«3 

XM 

KM 

FRAC 

FRAC 

32. 

A 

13.332 

0. 

192 

12.810 

6.048 

6.233 

0*f41 

0.ft7 

32. 

2 

13.400 

0. 

192 

13.228 

6.042 

6.221 

0.013 

0.026 

32. 

1 

13.876 

0. 

195 

13.443 

6.039 

6.214 

0*032 

0.028 

32. 

0 

13.995 

0. 

219 

13.661 

6.037 

6.208 

0.024 

0.030 

31. 

6 

14.413 

0. 

195 

14.572 

6.025 

6.184 

-0.011 

0.027 

31. 

5 

15.439 

0. 

206 

14.810 

6.022 

6.178 

0.043 

0.027 

31. 

3 

14.947 

0. 

199 

15.297 

6.017 

6*165 

-0.023 

0.029 

31. 

2 

15.153 

0. 

203 

15.548 

6*014 

6.159 

-0.025 

0.027 

31. 

0 

15.627 

0. 

208 

16.061 

6.008 

6.147 

-0.027 

0.028 

30. 

8 

16.701 

0. 

221 

16.593 

6.002 

6.135 

0.007 

0.027 

30. 

7 

16.723 

0. 

221 

16.866 

6.000 

6.129 

-0.008 

0.020 

30. 

5 

17.451 

0. 

227 

17.426 

5.994 

6.117 

0.001 

0.017 

30. 

4 

17.518 

0. 

255 

17.713 

5*991 

6.112 

-0.01  1 

0.014 

30. 

0 

19.096 

0. 

243 

18.914 

5.980 

6.088 

0.010 

0.008 

29. 

9 

19.421 

0. 

253 

19.227 

5.977 

6.082 

0.010 

0.009 

29. 

7 

19.072 

0. 

243 

19.870 

5.972 

6.071 

-0.040 

0.020 

29. 

6 

19.144 

0. 

245 

20.200 

5.969 

6.065 

-0.052 

0.031 

29. 

5 

19.605 

0. 

258 

20.536 

5.966 

6.059 

-0*045 

0.036 

29. 

3 

20.929 

0. 

296 

21 .226 

5.960 

6.047 

-0.014 

0.037 

29. 

1 

20.666 

0. 

283 

21.941 

5.955 

6.036 

-0.058 

0.045 

29. 

0 

20.896 

0. 

294 

22.307 

5.952 

6.030 

-0*063 

0.050 

28. 

9 

23.378 

0. 

444 

22.681 

5.949 

6*025 

0*031 

0.046 

28. 

6 

21.589 

0. 

327 

23.840 

5.941 

6.008 

-0.094 

0.059 

28. 

5 

23.648 

0. 

400 

24.241 

5.938 

6.002 

-0.024 

0.060 

28. 

3 

24.303 

0. 

439 

25.063 

5.933 

5.991 

-0.030 

0.055 

28. 

2 

25.138 

0. 

437 

25.485 

5.930 

5.985 

-0.014 

0.048 

28. 

1 

26.962 

0* 

433 

25.914 

5.927 

5.980 

0.040 

0.050 

27. 

9 

27.242 

0. 

445 

26.797 

5.922 

5.969 

0.017 

0.027 

27. 

7 

28.263 

0* 

498 

27.711 

5.917 

5.957 

0.020 

0.026 

27. 

6 

27.341 

0. 

449 

28.180 

5.914 

5.952 

-0.030 

0.026 

27. 

4 

30.125 

0. 

599 

29.144 

5*909 

5.941 

0.034 

0.029 

27. 

2 

32.794 

0. 

599 

30.143 

5.903 

5.930 

0.088 

0.046 

27. 

1 

30.542 

0. 

609 

30.655 

5.900 

5.925 

-0.004 

0.045 

27. 

0 

32.576 

0. 

598 

31.178 

5.898 

5.919 

0.045 

0.049 

26. 

9 

31.091 

0. 

558 

31.709 

5.895 

5.914 

-0.019 

0*047 

26. 

7 

35.629 

0. 

632 

32.801 

5*890 

5.903 

0.086 

0.059 

26. 

6 

35.194 

0. 

620 

33.361 

5.887 

5.898 

0.055 

0.051 

26. 

5 

32.020 

0. 

588 

33.932 

5.884 

5.892 

-0.056 

0.057 

26. 

3 

37.662 

0. 

671 

35.105 

5.879 

5.882 

0.073 

0.062 

26. 

2 

37.429 

0. 

669 

35.707 

5.876 

5.876 

0.048 

0.065 

TOT.  RMS  FRAC.  RHOMEAS-RHOCALC  DIFP-  0.0AA75 
TOT.  RMS  FRAC.  RHOMEAS.  UNCERTAINTY^  0.03909 
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Density  (g/m 


Density  (g/nn 


Examination  of  the  variation  of  RMS  FRAC  with  altitude  shows  those 
regions  in  which  the  best  fit  is  obtained.  At  the  highest  altitudes  the  aver- 
age difference  for  any  single  point  is  as  high  as  7-8%,  due  to  the  small 
nxunber  of  counts  recorded  there.  As  the  sonde  descends  it  goes  through 
a region  where  the  average  difference  is  on  the  order  of  2%.  There  is, 
however,  a lower  altitude  region  (below  about  20  km)  in  which  beta  par- 
ticle energy  loss  dominates  the  effects  of  scattering,  and  the  count  rate 
decreases  as  z decreases  (Ref.  3.3).  Thus,  near  20  km  there  is  little 
variation  of  coimt  rate  with  z,  and  the  density  uncertainty  becomes  very 
large  in  spite  of  the  large  count  rates.  This  region  is  being  approached 
at  the  lowest  altitudes  in  Table  3.5,  which  accounts  for  much  of  the  aver- 
age difference  of  about  5%  there.  As  noted  previously,  the  4 measure- 
ments made  in  the  vicinity  of  the  64  km  apogee  during  a period  of  about 
20  seconds  are  particularly  important,  since  the  sonde  descends  rapidly 
for  some  distance  thereafter.  The  data  from  these  four  points  could  have 
been  averaged  together  to  give  one  measurement  at  the  average  altitude  of 
64.0  km.  The  re.'ivilt  wovild  have  been  0.1833  g/m^  with  a fractional  uncer- 
tainty of  .054  (about  half  that  of  each  of  the  four  points).  The  calculated 
result  from  the  fit  is  .190  g/m^,  a fractional  difference  of  only  .035. 

Thus,  these  four  points,  as  utilized  with  least  square  procedure,  serve 
to  define  p(z)  accurately  in  the  apogee  altitude  region. 

It  is  of  interest  to  compare  the  H(z)  values  with  those  of  an  approp- 
riate model  atmosphere,  such  as  that  of  the  Spring/Fall  Mid- Latitude  Model 
(Ref.  3.1).  Values  of  H(z)  are  not  given  in  that  particular  model.  It  is 
straightforward,  however,  to  derive  the  associated  H(z)  by  fitting  the  two 
p(z)  values  that  bracket  the  altitude  of  interest.  By  use  of  this  procedure 
and  the  results  of  Table  3.5  we  obtain  the  following: 


Table  3.6 

Comparison  of  Scale  Heights 
for  4/23/76  Flight 


“model"'"''  “fit"'"'* 


25 

6.35 

~5.  85 

30 

6.  50 

6.  09 

35 

6.  47 

6.40 

40 

6.85 

6.  76 

45 

7.?4 

7.  18 

50 

8.04 

7.67 

55 

8.40 

8 24 

60 

8.  11 

8.  94 

65 

8.08 

~9.  70 

25 


In  general  the  H(z)  values  from  the  analytical  fit  are  slightly  below  those 
of  the  model  in  the  25-30  km  region,  close  in  the  35-55  km  region,  and 
somewhat  higher  in  the  60-65  km  region.  Certainly,  the  agreement  is 
very  reasonable  on  the  average. 


3.5  Summary  and  Discussion 


In  summary,  we  believe  that  the  results  of  applying  the  least  square 
mathematical  procedure  developed  in  Section  3.2  to  the  LBB  flight  data  in 
Section  3.3  show  the  procedure  to  be  valid.  An  accurate  analytical  fit  was 
obtained  that  is  capable  of  representing  the  entire  tabular  array  of  data, 
and  we  believe  that  the  procedure  should  be  applied  to  future  LBB  or  Beta- 
sonde  II  flights. 


It  should  be  observed  that  a similar  procedure  can  also  be  applied  to 
other  altitude-dependent  measurements,  for  example  temperature  and  pres- 
sure, provided  a tabular  array  of  the  measured  quantity,  and  its  uncertainty, 
is  available.  In  the  case  of  pressure  measurements  the  above  procedure, 
could  be  used  almost  in  its  exact  form,  since  pressure  also  depends  expo- 
nentially on  altitude.  For  temperature  the  analytical  function  would  have 
to  be  chosen  to  represent,  in  an  average  manner,  the  known  variation  of 
that  quantity  with  altitude. 


4.  Atmospheric  Density  Response 


In  this  section  the  Betasonde  atmospheric  density  response  is  formu- 
lated from  an  analytical  viewpoint,  and  the  relative  responses  of  various 
sources  are  derived  by  use  of  the  formulation.  This  then  allows  a choice 
of  the  optimum  source  for  the  sonde,  and  a calculation  of  the  expected  accuracy. 

4.1  Form\xlation 


Consider  the  beta  particle  forward  scattering  configuration  shown  in 
Fig.  4.1. 


Collision  with 
atmospheric  atom 


Emitted 
beta  particle 


\5X^ 


Annular  Source 


Shield 


Scattering 

angle 


"W 

Detector 


Fig.  4.1  Beta-ray  Scattering  Configuration 


The  single- scattering  count  rate  Cg(E^.)  cps  of  the  detector  for  the  source 
s and  detector  energy  threshold  E^.  (the  energy  above  which  all  electrons 
are  counted)  is  found  by  integrating  over  all  scattering  angles  0 into  solid 
angles  dQ: 

00 

Cg(E^)  = p (nN^/M)  J f S^(E^)Bg(E)da(E,e)dE  (4.1 


Here  p = atmospheric  density,  g/cm^ 

= Avogadro's  number,  6.023x10^^  molecules/(g-mole) 
M = molecular  weight  of  atmosphere,  g-mole 
n = number  of  atmospheric  atoms /molecule 
da  = f(0)F(E)dJ2  cm^/atom 
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iW'j' I 'i* 


F(E)  = (l-p2)/p4 

f(0)  = (Z2ro2/4)  Sin-40/2  cm^/(sr-atom) 

P = v/c  = ratio  of  beta  particle  velocity  to  that  of  light 
Z = atomic  number  of  atmosphere 
rQ  = classical  electron  radius,  2.82xl0“^^cm 
Sg  (E^)  = total  emission  rate  of  source  for  particles  of  energy 
greater  than  E^.*  P/sec 

Bg(E)  = spectrum  for  beta  source  of  type  s,  normalized  to  unity 
from  Ej.  to  00,  p/keV 
E = kinetic  energy,  keV 

The  cross  section  dcr  (E,  0)  for  scattering  of  an  electron  of  energy  E (in  the 
energy  range  of  interest  here,  550  keV)  through  an  angle  0 is  given  by  Evans 
(Ref.  4.1).  As  shown  above,  it  is  separable  into  two  functions,  the  first  de- 
pending only  on  E,  and  the  second  only  on  0.  The  function  F(E)  can  be 
written  conveniently  in  terms  of  the  energy  E as 


F(E) 


1 + R 


(1  + R)-^  - 1 


where 


R = E/(m^c2) 


and  =511  keV  ( 

is  the  electron's  rest  mass  energy. 

Because  the  energy  and  angular  integrations  can  be  separated,  we 
write  the  count  rate  in  terms  of  the  single  scattering  constant  Kg, 

Cg(E^)  = K^(E^)p  cps,  ( 


where 


K (E.)  = GS  (E.)  e (E.)  (cps-cm^)/g, 

S V W Ir  S V 


(4.6) 


G is  a constant  for  a given  geometry  and  gas, 

G = {nNjU)J'm)dil 
« f2 

and  the  "beta  source  efficiency  factor"  is 


g-P  ’ 


(4.7) 


OO 

^s^Et)  = f (E)F(E)dE. 
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I 

I 

For  any  particular  beta  source,  the  spectrum  Bg(E)  is  either  known  or  can  ; 

be  measured.  Hence  the  efficiency  factor  £g{E|.)  can  be  readily  determined  i 

for  any  threshold  energy  E^.  Suppose  for  that  particular  source  the  constant  i 

Kg{E^)  has  been  determined  by  measuring  the  slope  of  the  count  rate  versus 
density  curve  for  the  gas  of  interest.  Then  the  constant  G for  that  geometry 
and  gas  can  be  calculated  from  | 

G = Kg(E^)/[S3(E^)£g(E^)]  (4.9) 

As  a consequence,  for  any  other  source  the  value  of  Kg(E^)  can  be  deter- 
mined for  that  gas  and  any  desired  energy  threshold  from  (4.6). 

4.2  Beta  Source  Evaluation 

From  (4.5)  it  is  clear  that  the  larger  is  Kg(E^)  the  higher  will  be  the 
count  rate  for  any  density  p.  Kg(E^)  depends  on  E^.  through  both  the  efficiency 
factor  and  the  emission  rate  Sg(E|.)  of  betas  having  energy  greater  than  Et* 

It  is  convenient  to  write  this  as 

Ss(E^)  = fs(E^)S^o 

I 

I 

where  S^q  is  the  total  emission  rate  (p/sec)  and  fg(E^.)  is  the  fraction  of  those  | 

betas  that  have  energy  greater  than  E.  for  the  source  type  s.  I 

i 

In  the  present  instrument  the  energy  threshold  E^  must  be  set  at  about  | 

50  keV  in  order  to  exceed  the  electronic  noise  sufficiently.  This  means  that  | 

a beta  source  such  as  Ni-63,  with  a maximum  energy  of  67  keV,  will  not  | 

have  a large  value  of  f(E^).  However,  the  fact  that  the  maximum  energy  i£^  | 

only  67  keV  causes  £(E^)  to  be  large.  Thus  we  rewrite  Kg(E^)  as  follows:  S 


K (E^)  = GS  Q (E J 
s'  t so  s'  t 

1 

(4.11) 

where 

(4.12) 

measures  the  overall  "quality"  of  the  source  of  type  s for  the  energy  thresh- 
old Ej..  Thus 

Q^(E^)  = C(E^)/(GpSgo)  (4.13) 

determines  the  count  rate  per  unit  of  total  source  activity  and  per  unit  den- 
sity for  a fixed  energy  threshold  E^  and  geometry. 
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By  use  of  a tabulation  of  beta  spectra  (Ref.  4.Z)  for  Tl-204,  Sr-Y-90, 
and  Ni-63,  and  an  experimental  measurement  for  Pm-147,  the  results  given 
in  Table  4.1  were  obtained  for  = 50  keV.  This  is  a tabulation  of  Bg{E), 
F(E),  and  £{E)  = Bs(E)F(E)  versus  energy  for  each  of  these  sources.  At 
the  end  of  each  tabulation  is  the  result  obtained  for  from  (4.8).  By 

use  of  these  results,  and  the  fractions  fs(Ej.)  found  from  Ref.  4.2,  we  then 
obtain  the  results  in  Table  4.2. 
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Table  4.2 

Atmospheric  Density  Measurement 
Evaluation  Factors,  E^  = 50  keV 


Source  Type 

Half-life  (yrs) 

fs(E^) 

£s(Et) 

Qs(E^) 

Thallium -2  04 

3.6 

. 14 

3.  24 

.45 

St  r ont  ium  - Ytt  r ium  - 9 0 

28 

.95 

2.  27 

2.  16 

Nickel-63 

120 

. 02 

24.45 

.49 

Pr  om  ethium  - 1 4 7 

2.  5 

. 50 

12.  09 

6.  05 

Ihe  higher  is  the  value  of  Qg  the  smaller  is  the  amount  of  source  activity  | 

required  to  obtain  the  same  count  rate  (i.  e.,  density  accuracy)  at  a given  j 

density.  From  that  viewpoint,  Pm-147  is  the  best  choice.  Unfortunately,  | 

it  has  the  shortest  half-life  and  would  probably  require  calibration  of  a I 

given  sonde  at  least  on  an  annual  basis.  It  is,  however,  available  at  a I 

reasonable  price  (less  than  about  $1000)  in  higher  activity  amounts  (up  to  j 

about  a curie  in  the  annular  configuration  used  here)  than  any  of  the  others.  j 

If  it  were  possible  to  develop  a detector  with  a much  lower  energy  thresh- 
old (;$,10-20  keV),  then  the  Ni-63  source  would  definitely  be  optimum.  In 
that  case  fg(Et^  would  be  near  0.50,  and  £s(E|.)  would  be  even  larger. 

This  would  be  necessary,  because  the  source  is  not  readily  available  in 
amounts  larger  than  50-100  mCi. 

For  the  present  system  we  believe  that  Pm-147  is  optimum.  It  will 
provide  a means  of  establishing  the  utility  of  the  betasonde  for  the  upper 
portion  of  the  middle  atmosphere. 

4.3  Estimated  Density  Accuracy 

The  source  strength  of  the  same  annular  Pm-147  beta  source  pre- 
viously flown  on  the  Low  Background  Betasonde  (LBB)  has  been  measured  j 

on  May  12,  1 977  to  be  1 

Sg(E^)  = 2.  76  X 10*^  p/sec  = 74.5  mCi  (4.14) 

' (test  source) 

for  E^  = 50  keV.  The  slope  of  the  count  rate  versus  density,  Eq.  (4.5),  was 
measured  in  the  geometry  of  Fig.  2.*1  to  obtain  the  calibration  constant 

\ I 

9 3 ; 

Kg(Ep  = 0.93  X 10^  cps/(g/cm^)  (4.15)  ; 

(test  source)  [ 

for  the  source  intensity  (4.14).  ! 
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A difference  between  the  present  sonde  and  the  LBB  is  that  a signif- 
icant wall  effect  is  present  during  calibration  in  the  Panametrics  chamber 
(about  2'x3').  This  was  not  experienced  by  the  LBB  due,  we  believe,  to 
the  extremely  collimated  geometry  of  the  geiger  counters  themselves. 

Thus,  going  to  the  semiconductor  detectors  causes  a significant  increase 
in  overall  efficiency,  but  also  produces  a wall  effect  in  the  small  chamber. 
This  does  not  affect  the  determination  of  the  constant  Kg.  since  it  is  the 
slope  of  the  count  rate  versus  density  curve.  It  does  mean,  however,  that 
in  order  to  obtain  an  actual  calibration  curve  at  densities  less  than  that 
equivalent  to  about  70  km,  it  will  be  necessary  to  use  a much  larger  cham- 
ber. Such  chambers  are  available  at  NASA  Langley  Research  Center  (both 
20'  and  60'  diameter  spheres)  and  have  been  used  previously  (Ref.  3.3)  for 
betasonde  calibration  work. 

Once  it  is  established  that  a calibration  measurement  Kg  in  a small 
chamber  leads  to  the  same  result  as  that  obtained  in  a large  chamber  - in 
which  the  wall  effect  is  much  reduced  - it  is  not  then  necessary  to  recalib- 
rate in  the  large  chamber  unless  the  source-detector  geometry  is  changed 
in  some  significant  way.  Rather,  a value  of  Kg  can  be  carefully  measured 
in  the  small  chamber  in  the  z <60  km  region,  where  the  wall  effect  is  only 
a fraction  of  total  count  rate,  even  in  the  small  chamber,  and  the  calibration 
curve  for  higher  altitudes  can  be  obtained  by  extrapolation.  This  is  possible 
because  single -scattering  is  by  far  the  dominant  effect  for  z^40  km,  hence 
the  calibration  curve  is  extremely  linear  and  the  count  rale  versus  density 
is  given  rigorously  by  Eq.  (4.5). 

Thus,  the  calibration  results  in  the  small  chamber  can  be  used  to 
estimate  the  accuracy  obtainable  at  high  altitudes  and  with  a more  intense 
source. 


The  fractional  statistical  uncertainty  for  any  count  rate  Cg  measured 
for  a time  t(sec)  is  given  by  (Ref.  3.2) 


(4.16) 


Ap/p  = i/yct 


By  use  of  (4.5)  this  is 


Ap/p  = i/yK^(E^)p, 


(4.17) 


This  shows  explicity  that  the  statistical  uncertainty  is  a function  of  the 
threshold  Ej.  through  the  calibration  constant  Ks(Et).  Of  course.  Kg  is  also 
proportional  to  the  source  strength  Ss(Ej)  through  (4.6).  We  believe  that  a 
source  of  intensity  about  5 times  higher  than  (4.14)  should  be  obtained  to  use 
as  the  first  flight  source, 
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for  which 


S (E^)  = 1.38  X 10^°  p/sec  = 373  mCi, 

(flight  source) 


K^(E^)  = 4.65  X 10^  cps/(g/cm^). 

(flight  source) 


(4.18) 


(4.19) 


Now,  by  use  of  the  standard  atmosphere  densities  (Ref.  4.3),  Eqs-  (4.5), 
(4.19)  and  (4.17),  the  following  results  are  obtained: 


Table  4.3 

Estimated  Fractional  Statistical  Uncertainty  for  Betasorde  II 
and  First  Flight  Source  (S  « 373  mCi) 


z 

km 

^ 3 

g/cm-’ 

Cs 

cps 

Ap/p 

t = 20  sec 

t = 5 sec 

t = 1 sec 

60 

3.1  xlO-"^ 

1441 

.006 

.012 

.026 

70 

8.8x10-8 

409.2 

.011 

.022 

.049 

80 

2.0x10-8 

93.0 

.023 

.046 

.104 

90 

3.2x10-9 

14.9 

.058 

.116 

.259 

100 

5.0x10-10 

2.3 

.147 

.295 

.659 

The  uncertainties  listed  in  this  table  are,  of  course,  the  minimum  that  will 
occur  for  the  indicated  coxint  rates.  The  actual  total  error  will  be  slightly 
larger  due  to  some  (as  yet  undetermined)  cosmic  ray  effect.  This  is  ex- 
pected to  be  extremely  small,  however;  much  less  than  in  the  LBB  (Ref. 

3.2). 

As  seen  from  Table  3.1,  the  betasonde  spends  about  20  seconds  with- 
in a 1 km  region  near  apogee.  As  discussed  above,  this  allows  a very  ac- 
curate measurement  of  the  density  in  this  particular  region;  the  procedure 
described  in  Section  3 provides  an  analytical  density  fit  that  takes  this  into 
account.  Th-  RMS  deviation  of  the  fit  from  the  measured  results  approaches 
the  statistical  accuracy  both  in  that  high  altitude  region  and  on  an  overall 
average  basis  (Table  3.5).  Thus,  it  can  be  expected  that  the  accuracy  of 
the  lit  would  approach  about  +6%  near  90  km  for  an  apogee  in  that  region. 

At  lower  altitudes,  as  the  descent  velocity  increases,  it  was  found  that  5 
second  intervals  provided  a meaningful  measurement.  Hence  from  about 
80  km  dovnward  to  perhaps  40-50  km  the  accuracy  will  be  no  worse  than 
about  +5%,  as  seen  in  Table  4.3. 
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5.  Phase  I Summary 

During  Phase  I,  Research  and  Development,  the  following  work  has 
been  accomplished. 

1.  A technique  has  been  developed  for  making  an  analytic  fit  to  an 
entire  array  of  density  versus  altitude  data  on  a least- square- 
basis.  The  procedure  was  applied  to  data  from  a Low  Back- 
ground Betasonde  flight  with  good  results. 


2.  By  use  of  theoretical  calculations  and  experimental  measure- 
ments, it  was  shown  that  Pm-147  is  the  optimum  choice  source 
for  the  single- scattering  betasonde.  This  conclusion  would  be 
modified  if  it  were  possible  to  reduce  the  energy  threshold  (that 
energy  above  which  all  betas  are  counted),  to  about  10-20  keV. 
In  that  case  Ni-63  would  probably  be  optimum. 

3.  A laboratory  model  of  Betasonde  II,  the  semiconductor -detector 
version,  has  been  constructed  and  calibrated  with  Pm- 147.  Al- 
though some  wall- effect  count  rate  was  evident  at  very  low  den- 
sities (due,  apparently,  to  the  open  geometry  of  the  detector), 

it  was  possible  to  obtain  meaningfiil  calibration  data. 

4.  The  calibration  data  were  used  to  show  that  with  a Pm-147 
source  of  about  500  mCi,  it  should  be  possible  to  obtain  about 
+6%  accuracy  in  the  90  km  region  for  an  Areas  launch  having 
apogee  near  that  altitude.  Below  that  altitude  the  accuracy 
would  be  better. 

5.  The  flight  instrument  must,  ultimately,  be  calibrated  in  a large 
{40'x60')  chamber  in  order  to  verify  the  high  altitude  (>70  km) 
portion  of  the  calibration  curve. 
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During  Phase  II,  Design  and  Fabrication,  the  flight  unit  v.'ill  be  com- 
pleted and  tested.  It  will  be  calibrated  by  use  ot  the  nev.  ~500  mCi  Pm-147 
source,  and  the  unit  will  be  ready  for  delivery  at  the  completion  of  the  work. 
During  the  final  phase  of  the  work  (not  presently  expected  to  be  funded  by  the 
Army  Research  Office)  the  instrument  rn<..3t  be  integrated  into  an  Ar^as  pay- 
load  for  launch  at  White  Sands  Missile  Range.  The  Betasonde  II  should  then 
be  flown  to  the  80-90  km  region  at  least  twice,  and  the  analytical  fitting  pro- 
cedure presented  above  shoxild  be  applied  to  the  data.  This  will  form  the 
basis  for  routine  direct  measurement  of  atmospheric  density  up  to  at  least 
the  80-90  km  region.  Should  this  prt>ve  to  be  feasible,  as  is  expected  based 
on  the  present  results,  it  would  then  be  useful  to  consider  application  of  the 
technique  up  to  the  100  km  region  either  by  use  of  a more  intense  source 
(several  kCi),  by  developing  detection  techniques  allowing  a lower  energy 
threshold  (<10-20  keV),  or  a combination  of  these  measures. 
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